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The bilayer quantum Hall system at a total fiUing of i/t = 1 has long resisted explanation in terms 
of a true counterflow superfluid, though many experimental features can be seen to be "almost" that 
of a superfluid. It is widely believed that quenched disorder is the root cause of this puzzle. Here we 
model the nonperturbative effects of disorder by investigating the v — 1 bilayer in a strong periodic 
potential. Our model assumes that fermions are gapped and real spins are fully polarized, and 
concentrates on the pseudospin variable (the layer index), with the external potential coupling to the 
topological (Pontryagin) density of the pseudospin. We find that as the potential strength increases, 
there are ground state transitions in which the topological content of the pseudospin configuration 
changes. These transitions are generically weakly first-order, with a new quadratically dispersing 
mode (in addition to the linearly dispersing Goldstone mode) sometimes becoming nearly gapless 
near the transition. We show that this leads to strong suppressions of both the Kosterlitz-Thouless 
transition temperature and the interlayer tunneling strength, which we treat perturbatively. We 
discuss how these results might extend to the case of true disorder. 

PACS numbers: 73.50.Jt 



I. INTRODUCTION 

It has long been expected theoretically that bilayers 
offer thepossibility of Bose condensation of interlayer 
excitons[l[. While investigations in electron- hole bilay- 
ers are ongoing Q, the electron bilayer system in a large 
magnetic field at total electron filling ~ 1 offers a differ- 
ent way of realizing this goalQ . There are two important 
length scales in this system, the interlayer separation d 

and the magnetic length / = y^^, (MKS units) where e 
is the magnitude of the electron charge and B is the mag- 
netic field. In addition to the dimensionless tuning pa- 
rameter d/l, there may also be a layer imbalance v-^ — v^, 
where f: i refer to the top and bottom layer, and an 
in-plane field . The discussion in this paper will be re- 
stricted to the case of balanced layers with 1/^ = 1^1 = ^ 
and = 0. In this case, by a particle-hole transforma- 
tion in one layer, one can immediately see the possibility 
of exciton formation and condensation for small d/l. 

The first theoretical works on this system two decades 
ago predicted in the clean limit that at T = and small 
d/l, the system could be thought of either as a exci- 
ton Bose condensate or as a pseudospin quantum Hall 
ferromagnet U with a planar anisotropy. The presence 
of a gapless, linearly dispersing, Goldstone mode (in the 
absence of interlayer tunneling) was notedQ indicating 
a spontaneous broken symmetry, as was the expectation 
that when the interlayer tunneling amplitude h became 
nonzero but small compared to the other energy scales, 
there would be a delta-function peak in the interlayer 
conductance Gjl at zero biasQ. Upon the application 
of an in-planc field , the peak is expected to split into 
two symmetric peaks separated by a bias voltage propor- 
tional to 0. Furthermore, the tunneling is expected to 



occur within a Joscphson length {Ij 



with J being 



the pseudospin stiffness) of the contacts@,01- Finally, for 
h — 0, a, Bcrczinskii-Kostcrlitz-Thouless transition[3, 0| 
at a nonzero temperature Tkt is expected, and the coun- 
terfiow conductance, where current fiows in opposite di- 
rections in the two layers, is expected to be infinite for 
T < Tkt^- An excellent review summarizes the theo- 
retical expectations in the clean svstem[loj. 

Experimentally, it is found that as d/l varied, there is 
a transition[ll| between a compressible phase with un- 
quantizcd Hall conductance for large d/l (presumably 
adiabatically connected to a system with two decoupled 

= ^ Composite Fermion[l^ Fermi seas[l3 in the two 
layers) and an incompressible phase for d/l < 1.6 with a 
quantized Hall conductance. 

Where the experimental results deviate from theoreti- 
cal expectations is in the interlayer tunneling [T^ - [l6| and 
counterflowfirj properties of the incompressible phase. 
It is never possible in an experimental sample to make h 
precisely zero, but it can be made much smaller (of order 
lOOfiK) than any other energy scale, including the tem- 
perature. In such samples, there is a peak at zero-bias 
in Gil, but its width remains nonzero even at the lowest 
temperatures [1^ . Upon the application of an in-plane 
field, features at nonzero bias are seen to vary in a linear 
manner with Sy, but the peak at zero-bias remains[l6[, 
contrary to theoretical expectations, and retains most of 
the weight. The counterfiow conductivity increases as T 
decreases, but remains finite down to the lowest measured 
T[17|. Contrary to theoretical expectations, interlayer 
tunneling occurs throughout the area of the sample [18 1 
rather than within a Joscphson length of the contacts. 
Recently, critical interlayer tunneling currents have been 
measured [l^ (with some puzzling aspects [l^l), and are 
also found to be proportional to the area of the sam- 
ple. Finally, while a sharp drop in the peak value of the 
interlayer conductance is seen at a particular T, the phe- 
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nomenology of this nonzero-temperature transition[2l[ (if 
indeed it is a thermodynamic transition) does not seem 
to be Bcrezinskii-Kosterhtz-Thouless-hke. 

Much theoretical work0-0: l22| - |43j has been carried 
out on this system in the past two decades in several 
directions. We will not attempt to review the vast litera- 
ture, except to comment that no satisfactory fundamen- 
tal explanation of the entire spectrum of discrepancies 
noted above has been proposed yet. It has been recog- 
nized that an extra mechanism of dissipation seems to 
be active in this system, as seen in the finite value of the 
counterflow conductivity [13], and in the nonzero width of 
the zero-bias peak in Gil- Assuming such a mechanism 
phenomenologically, several authors have successfully de- 
scribed some fraction of the experimental results, as for 
example in Ref. [Ill- Recent efforts have focused on 
the effects of quenched disorder [13 [H [H, js^, US [H- 

, and succeeded in explaining some aspects of the ex- 
periments. Theoretical evidence from studies of critical 
counterflow [111 and tunneling currents [i^ also lead to 
the identification of quenched disorder as the root cause 
of the discrepancies between theory and experiment. 

In this paper we will start with the assumption that 
quenched disorder is responsible for the qualitative differ- 
ences between theory in the clean limit and experiment. 
As Efros 4J] showed long ago, in any incompressible sys- 
tem, disorder necessarily has a nonperturbative effect, 
because the clean system cannot screen linearly. In two- 
dimensional electron gases, there is the further feature 
that the fluctuations of the disorder potential, created 
by positional flutuations of the remote dopants, diverge 
at long distances. In the semiclassical Efros picture, the 
quantum Hall system screens by forming compressible 
puddles (of typical size s, the distance between the 2DEG 
and the dopant layer) separated by incompressible strips 
of width a few magnetic lengths. This picture is sup- 
ported by imaging experiments [4^. Based on this pic- 
ture of strong smooth disorder, two of the present au- 
thors previously presented [s^ a classical coherence net- 
work model displaying some of the experimental phe- 
nomenology, such as the proportionality of the interlayer 
tunneling conductance to the areafisj rather than length 
of the contacts. A more detailed, though still classical, 
calculation [12 has recently used the coherence network 
model to predict a large enhancement of the Josephson 
length due to disorder. 

Here we take a slightly different approach, in- 
spired by studies of the Bose-Hubbard model[4^ and 
one-component quantum Hall systems in a periodic 
potential (47|. These studies were motivated by trying to 
develop a field-theoretic understanding of the superfluid- 
Mott Insulator transition for bosons or the plateau tran- 
sition for quantum Hall systems. Note that without a po- 
tential, Galilean invariance means that one cannot even 
obtain a transition in both these cases. Once a second- 
order transition has been obtained in the presence of 
a periodic potential, one may examine the relevance or 
irrelevance of operators pertaining to various kinds of 



disorder [48j at this fixed point. 

Combining ideas from the Efros picture [i^l and the 
plateau transition work[47j. one may hope that a peri- 
odic potential will capture the essential nonperturbative 
features of quenched disorder, and true disorder can be 
added later [48|. In this paper we will carry out the first 
part of this program, ending with some speculations re- 
garding true quenched randomness at the end. 

As a further simplification, we will assume that 
fcrmionic degrees of freedom are gapped out and focus 
exclusively on the pseudospin degrees of freedom as the 
only relevant ones. This also involves freezing out the 
real spin, which may not be quite correct for typical 
experimental samples (iol. [sof. However, recently it has 
been shown that the interlayer coherent phase is robust 
to large Zeeman energies, and survives full polarization 
of real spin[5l|. Our study will be of direct relevance to 
that system, and the qualitative results will likely hold 
for the low Zeeman energy case as well. From now on, 
we will refer to "spin" always meaning the pseudospin 
vector, represented by a unit vector n. An important in- 
gredient is the spin-charge relation [52],[5^, which holds in 
multicomponent systems in a quantum Hall phase[7l[l0|. 
In our context, this states that slow variations of spin 
contain Coulomb charge: 



(1) 



It is through this 5p that the spins couple to the external 
potential. For ease of computation we will put our spins 
on a square lattice, with the lattice spacing chosen to 
be of the order of the magnetic length /. In this case the 
above, continuum, expression has to be replaced with the 
spherical area 5J| subtended by the three non-coplanar 
spins ni, n2, na, 



SQi23 = tan — ■ ■ 

27r V 1 + ni • n2 -I- n2 • ns + ns • ni 



(2) 



The Hamiltonian of our model has the form 



J E K(r)n,(r') + n,(r)n,(r')) + ^E(".(r)) 

(rr'> r 

-/i E (r) - E fix, Y)SQa + Hu [5Q^] . (3) 

r □ 



Here J is the nearest-neighbor ferromagnetic coupling, 
and is related to the continuum spin stiffness, F is the 
interlayer charging energy, h is the interlayer tunneling 
amplitude, and Vb is the strength of the periodic potential 
with /(X, Y) being its specific functional form, defined 
on dual lattice sites R = {X,Y). The sum in the third 
term, Ej refers to a sum over plaquettcs of the lattice, 
□ 

and we have also introduced a local Hubbard interaction 
Hjj to model the energy cost (of order the interaction 
scale Ec — ^) of localizing a charge within a distance 
of /. We will specify the precise form of these last two 
terms in the next section. Most of our results are for 
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J = r, with our Hubbard U being about ten times J, for 
different values of Vq. We also vary the size of the unit 
cell between 12 x 12 and 16 x 16 lattice sites, with each 
representing an area of order 100 — BOOP « (s is the 
distance between the 2DEG and the dopant layer), con- 
sistent with values of the expected size of Efros puddles 
for relevant experimental samples. 

Our results, which we summarize here, and elaborate 
upon in later sections, show that there are some surpris- 
ing features present in this model. Consider first h = 0. 
We find a sequence of ground state transitions where the 
topological content of the spin configuration changes. For 
sufficiently large U, the first transition is second-order 
while the others are generically weakly first-order. There 
is always a linearly dispersing Goldstone mode (which we 
will henceforth call the G-mode) indicating the broken 
U{1) symmetry for h = 0. The next higher mode is usu- 
ally quadratically dispersing (we will call it the Q-mode), 
with a gap that sometimes becomes very small near the 
transitions. In the case when the transition is contin- 
uous, the gap of the Q-mode vanishes at the transition 
continuously. It is important to note that true disorder 
will generically convert first-order transitions to second- 
order [H^, and thus the Q-mode is expected to be truly 
gapless at the ground state transitions in the disordered 
model. Two distinct effects can be traced to the transi- 
tions and the properties of the Q-mode near them. At 
values of Vq near the transitions, there are configurations 
close in energy with a different topological density, which 
means that the system becomes highly polarizable. In an 
effective XY model (to be described in Section) this leads 
to a greatly reduced core energy for vortices, or a greatly 
increased fugacity, which strongly suppresses the KT- 
transition temperature. Secondly, when the Q-mode gap 
is small, thermal fluctuations of the Q-mode are present 
at realistic temperatures, and strongly suppress the in- 
terlayer tunneling amplitude h. Finally, the low-lying 
Q-excitations offer one possible microscopic mechanism 
for the pervasive and puzzling dissipation seen in exper- 
iments. These are the main results of our work. 

Let us remark briefly about the realistic case h ^ but 
smaller than any other energy scale in the problem, in- 
cluding the temperature. The ground state phase transi- 
tions arc controlled by the competition between the core 
energies of merons/antimerons and the external poten- 
tial. Both these are on the scale of Ec- Therefore, we 
expect a nonzero but tiny h <^ T, V, Ec will not affect 
the transitions in any qualitative way. Similarly, for the 
collective mode dispersions, the primary effect of a tiny 
h is to gap the G-mode at q = 0. We expect all other 
modes to be robust against the introduction of a small 
h. 

It is important to note that in the clean system, h 
is relevant below Tj^tQ, and the regimes of small and 
large h are adiabatically connected to each other, with no 
phase transition separating them. This is not true for the 
system with a periodic potential or strong, smooth dis- 
order. As h is increased while all other parameters (such 



as Ec and V) are kept flxed, the system will undergo 
ground state phase transitions in which the topological 
density reduces, presumably becoming the uniform fer- 
romagnetic state for very large h ^ V . This qualita- 
tive distinction between weak and strong tunneling in 
the presence of a periodic potential (and presumably a 
strong smooth disorder potential) is consistent with ex- 
periments, which see the almost-superfiuid phenomenol- 
ogy only for weak h. 

The plan of this paper is as follows: In Section II 
we will present the precise definition of our model and 
present sample ground-state configurations. Section 
III presents our method for computing the collective 
mode spectra of excitations above the ground-state, and 
discusses the effects of the potential strength V , the 
inter-layer tunneling /i, and the interaction strength U 
on this. This is followed by a discussion of how one may 
extract the spin stiffness of our model from the numerical 
results. Section IV discusses two non-trivial fiuctuation 
effects in this system, first how the quadratically dispers- 
ing collective mode can suppress inter-layer tunneling, 
and then the suppression of Tkt near quantum phase 
transitions due to vortex excitations. We end with some 
conclusions, caveats, and speculations on what our stud- 
ies tell us about the case of true disorder in Section V. 



II. MODEL AND GROUND-STATES 

The form of the Hamiltonian of our model, Eq. ([3]), 
is guided by the following considerations: (i) We ignore 
fermionic excitations at the very outset, since in any in- 
compressible state, the energy to create a fermion is of 
the order of the Coulomb scale Ec, which is much larger 
than the other relevant energy scales. Similarly we ig- 
nore the dynamics of the real spin, assuming that it is 
fully polarized 5l| . (ii) We put the system on a lattice 
for computational convenience [28j . As long as all the 



relevant length scales, such as the size of the puddles, 
are much larger than the chosen lattice scale, we expect 
that our results will be physically correct. To be specific, 
we will choose the lattice spacing a = 1\/2tt such that 
each elementary plaquette has the area of one cyclotron 
orbit, and at z^t = 1, has one electron, (iii) In the con- 
tinuum, when there is no external potential, the Hamil- 
tonian should reduce to the continuum formQ in the 
long-wavelength limit, (iv) The coupling of the external 
potential to the spin (i.e., layer) degrees of freedom must 
be through the lattice generalizationlH^ of the topolog- 
ical charge [H, [5^. (v) To account for short-distance ef- 
fects beyond the spin stiffness, we introduce a Hubbard 
on-plaquette interaction term, (vi) Finally, we ignore 
the long distance component of the Coulomb interaction 
between induced charges for computational convenience, 
assuming that the Hubbard term can capture all quali- 
tative effects. 

We arc now ready to make the form of the Hamiltonian 
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precise. We consider a square lattice of points labelled 
by r, with the corresponding dual lattice labelled by R, 
with the convention that the site R corresponding to r 
lies half a lattice unit to the right and above r. The ex- 
ternal potential naturally lives on the dual lattice. Each 
dual site R defines a plaquette, whose sites can be la- 
belled ri, r2, Ta, r4, starting from r = ri and going 
counterclockwise. The topological charge corresponding 
to a given triplet of spins was given above in Eq. ([2|). 
The charge in a plaquette is defined as 



(5Qn(R) = SQ12: 



23 



1134 



^6Q 



124 



SQ 



234- 



(4) 



To fully specify the third term on the RHS of Eq. ^ 
we need the functional form of f{X,Y). We choose the 
simplest possible periodic form: 



= sinf-r^Uinf 



Na 



Na 



(5) 



Na is the size of the unit cell, which contains four pud- 
dles. We have experimented with other forms and ob- 
served no qualitative differences. Finally, we turn to the 
last term of Eq. ([3]) . In order to eliminate spurious con- 
figurations with equal and opposite topological charges 
in triplets 123 and 134 but no net charge in the plaque- 
tte, we write the Hubbard term as the symmetrized sum 
of all the triplets in a plaquette: 



= ( ('^Ql23)' + (5Ql34)' + (5Ql24)' + (<SQ 



R 



234) 



(6) 

Hartree-Fock calculations are a reasonable starting 
point for estimating the values of the various parameters 
entering the Hamiltonian. J and V are of the order of a 
few percent of Ec in such calculations, while U should be 
of the order of being the energy to localize a charge of 
order e within a length scale I. We will typically assume 
J = r, and U between 8 and 20 times J. Most of our 
results are for the size of the unit cell being iV = 16, but 
we will show results for TV = 12 as well. These values are 
realistic given I k, 200A and the distance to the dopant 
layer s sa 2000l. 

We parameterize the spins by z{y) = nz{r) 
and (j){r) which is the XY angle [ nx{r) — 
y/l — z{r)'^ cos0(r), ny{r) = \/l — z{r)'^ sin0(r) ]. We 
find the ground state configurations by starting with a 
random seed configuration, followed by simulated anneal- 
ing. Several different random seeds were tried at every 
value of Vo to eliminate the possibility of settling into a 
metastable minimum. 

We begin by showing the ground state spin configu- 
ration for a 16 X 16 unit cell at V = 3 and C/ = 8 (in 
units of J = r = 1). As shown in Fig. ([1]), there is a 
vortex/antivortex at the center of each puddle, partially 
screening the external potential. The length of the arrow 
indicates the projection of the spin in the xy-plane, while 
the color indicates whether the spin points in the up or 
down direction in z. A different way to plot the same 
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FIG. 1: {Color online) The ground state configuration for a 
16 X 16 unit cell with the strength of the periodic potential 
(in units where J — 1) being V = 3.0, and the Hubbard 
interaction is U = 8.0. The lengths of the spins denote their 
planar projection. Note a vortex/antivortex at the center of 
each puddle. 




FIG. 2: {Color online)The topological density of the spins in 
the ground state configuration for a 16 x 16 unit cell with 
V — 3.0, and U = 8.0 shown as a three-dimensional plot. 



configuration is to look instead at the topological den- 
sity of the spins. This is shown as a three-dimensional 
plot in Fig. ([2|). A more complicated configuration is 
shown in Fig. ([3]), where V has now been increased to 7. 
Each puddle now has two vortex-antivortex pairs. The 
corresponding topological density is plotted as a coun- 
tour plot in Fig. (|4]). The phase transitions correspond 
to values of V where the ground state energy of one con- 
figuration becomes lower than another. The very first of 
these transitions (for sufficiently large U) is second-order. 
For small V the ground state is uniform and ferromag- 
netic as at I^ = 0. At a critical V = Vc (which happens 
to be Vc = 4.2 for U — 18) there is an instability of the 
ground state, which, when followed, leads to a ground 
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FIG. 3: [Color online) The ground state configuration for 
V = 7.0 and U = 8.0. 




FIG. 4: [Color online) The topological cliarge in tlie ground 
state configuration for V = 7.0 and U = 8.0. 

state with a nonzero topological density. While the very 
first of these transitions is second-order, the rest seem 
to generically be weakly first-order. The ground state 
energies cross with different slopes as functions of V . 

In Fig. ([5]) we show the ground state of the 12x12 unit 
cell at = 3 and J7 = 8. The presence of two vortex- 
antivortex pairs is evident, as is the qualitative similarity 
to the case of the 16x16 unit cell of Fig. ([T]). Generally 
we find that all the qualitative features of the different 
ground states are the same for different sizes of the unit 
cell, but the specific values at which the transitions hap- 
pen do depend on unit cell size. 

III. DEVIATIONS FROM THE GROUND STATE 

Having obtained classical groundstates of this system, 
we are in a position to consider the energetics and dynam- 
ics associated with deviations from it. In what follows we 
will focus on small deviations, such that the system can 
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FIG. 5: [Color online) The configuration of the 12 x 12 lattice 
at V = 3 and U = 8. Comparing to Fig. ^ we see that the 
presence of the vortex-antivortex pairs is independent of the 
size of the unit cell, though the values of critical V at which 
the ground state transitions take place are different. 

be analyzed within a quadratic theory. Two types of 
quantities are of particular interest, the collective modes 
of the system, and the effective spin stiffness. The first of 
these can be used to explore the effects of thermal fluc- 
tuations, and the second is important in determining the 
effects of thermally generated vortices [i.e., Kosterlitz- 
Thouless physics.) 

A. Collective Modes 

Once a periodic ground state configuration z[r), 0(r) 
has been found, one can look at collective mode dynam- 
ics around this state. For dynamics, one needs the full 
action, which is 

S = Sb~ J dm, (7) 

where 5*^, the Berry phase term[55|. measures the spher- 
ical area covered by a closed path followed by a spin in 
time, 

1 

SB^-JdtJ d^^n(r,i,C)-atn(r,i,0 xacn(r,i,0. 

(8) 

Here ^ is an auxiliary variable and n(r, t, ^) is chosen such 
that at n(r,t, 0) = and n(r,t, 1) = n[t,r) [22|. Also, 
we have made the h explicit, and assumed that the spin 
is ^. This is reasonable since the underlying fermionic 
system has one electron per cyclotron orbit; with our 
choice of lattice spacing this translates to one electron 
per site. 

The Berry phase term may be conveniently expressed 

as 

Sb^ f dt\ ^ (z(r)^(r) - 0(r)i(r)) . (9) 
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Following standard methods, and recalling that any vari- 
ation of n has to be perpendicular to n we obtain the 
equation of motion 



dt U = n X 



6n 



(10) 



The ground state configuration is static, and satisfies 
dH , dU 



= 0. 



(11) 



dz{r)"'^ d4>{r)' 
Expanding in small fluctuations around the ground state, 
z(r) = z(r) + C(r), 0(r) = 0(r) + ^(r), (12) 
we obtain the linearized equations of motion 



As can be seen, the Hessian matrix of second derivatives 
(computed at the ground state configuration) appears on 
the right hand side. 



dw(r, s)dw{r' , s') 



(14) 



where we have introduced indices s, s' taking the val- 
ues 1, 2 for C and ip respectively, thereby organizing 
C(r) = w(r, 1), ip{r) ~ w{r,2) into a two-component 
vector. It is important to note that for h = one can 
rotate all the (j) by the same amount producing no change 
in energy, which means that the Hessian matrix H has 
an eigenvector with a zero eigenvalue. 

The next step is to exploit the larger unit cell and 
partially diagonalize the system in terms of a wave- vector 
k lying in the Brillouin Zone (BZ), arranging the 2N'^ 
variables in a unit cell into a single vector. We can now 
define the 2iV^ x 2iV^ matrices 



i^,,,s'(k:r,r') 



e-^^-i'^-'^')Hs.Ar,r'), (15) 

((^2) ss' ® Sr.r', (16) 



where (T2 is a Pauli matrix. It is henceforth understood 
that in all matrices r, r' are restricted to a single unit 
cell. Note that both matrices are Hermitian, but differ 
under transposition: 



{H{k)y = H{-k) 

= -K. 



(17) 



The earlier condition on the Hessian having an eigen- 
vector with zero eigenvalue translates to the vanishing 
of the smallest eigenvalue of as k — >■ 0. This will be 
important for the computation of the spin stiffness. 

Now the equations of motion can be written in a com- 
pact form 



wifV'-^'(k) =i?(k)v[-^i(k). 



(18) 



This is a generalized eigenvalue problem with both K and 
H being Hermitian. We first solve for the eigenvalues and 
eigenvectors of the auxiliary problem 



(19) 



to get the right eigenvectors . Taking the Hermitian 
adjoint and carrying out a few manipulations, one obtains 
the left eigenvectors 



t 



^r(k) X = C'(k). 



(20) 



Using the hermiticity of K and i/, and Eqs. ((T7)) and 
([20)) it is easy to show that the eigenvalues occur in pairs 
with equal magnitude and opposite sign. 

The normalization condition on the eigenvectors can 
be chosen to be 



(k) 



K^fik) 



sgn{ujx{k))5xx' , (21) 



where the sgn{u)\) is essential because K is not positive 
definite. Now, as a consequence of Eqs. ([TO]) and (PT|) 
the eigenvectors also satisfy 



^fl(k))'i7V'i^l(k) = K(k)|<5. 



AA' 



(22) 



Henceforth we will use only the right eigenvectors, drop- 
ping the superscript [R\. 

We find the following results for the collective mode 
dispersions as functions of V . Since we set the inter- 
layer tunneling /i = 0, there is always a gapless, linearly 
dispersing, Goldstone mode arising from the breaking of 
the continuous U{1) global symmetry. We will henceforth 
call this the G-mode. The G-modc is generically the low- 
est mode throughout the BZ, but may not be so close to 
a transition. The next higher energy mode is quadrati- 
cally dispersing at q = 0, and will play an important role 
in what follows. We will henceforth call it the Q-mode. 
At the first transition, which is second-order for large J7, 
the gap of this mode vanishes at q = 0. In Fig. ^ 
we show the dispersions of the two lowest energy modes 
at V — 4.2 and U — 24, close to the transition. It is 
clear that the Q-mode is the lowest energy mode for a 
substantial part of the BZ. 

At a generic transition, however, the gap of the Q- 
mode at q = does not vanish. Summarizing the be- 
havior of the Q-mode over the entire range of V we in- 
vestigated, we present Fig. ([7]), which shows the Q-mode 
gap as a function of for [/ ~ 8. We also examined 
the "wave functions" of the Q-mode near the transitions. 
The 0-component appears to be strongly peaked at the 
vortex/antivortex cores, while the z-component is more 
broadly distributed. Roughly speaking this may be un- 
derstood in terms of motion of the merons around their 
minimum potential sites, with a very flat effective curva- 
ture, indicating that "room" for further merons is de- 
veloping in the potential well. (Similar behavior has 
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FIG. 6: (Color online) The two lowest collective modes for 
V — 4.3 and U = 18. This is just before the transition 
from the uniform ferromagnetic state to a striped state, and 
is generically second-order at sufficiently large U. Note the 
almost gapless quadratically dispersing mode. 
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FIG. 8: The spinwave velocity of the G-mode as a function of 
V. Note that it never vanishes, but that it can vary by a factor 
of three, and have significant discontinuities at transitions. 

B. Ground-State Spin Stiffness 



V 



FIG. 7: The gap at q = of the Q-mode as a function of 
V. Note that it does not generically vanish at a transition, 
though there are transitions where it becomes particularly 
small. 



been observed in simulations of this system using an 
XY model [1^.) Moreover, the Q-modc dispersion is 
typically quite flat, indicating that these are fairly well- 
localized excitations. Another interesting point to note, 
shown in Fig. ([5]), is that there are large fluctuations 
of the G-mode velocity as V changes. It should be kept 
in mind that these are results to quadratic order in the 
fluctuations around the ground state configuration. In- 
teractions between the G- and Q- modes could affect the 
gap of the Q-mode and the spin- wave velocity as well. 



Another related quantity of interest for this system 
is the spin stiffness, which determines among other 
things the interactions between highly separated vortex- 
antivortex pairs. We begin by describing how it may be 
computed for an arbitrary ground state configuration. 
We define the spin stiffness by calculating the energy U 
for a particular spatially dependent configuration of (p 
and comparing it to the reference energy Ucont of the 
same configuration in a free continuum model, 



(23) 



The least expensive way to twist the system is controlled 
by the smallest eigenvalue of the Hessian H(k) as |k| — 
0. Note that these eigenvalues are not the same as the 
energies of the eigenmodes found in Eq. (jl9|) . which are 
the solutions to the dynamical problem. We consider a 
sample with lengths L^, Ly in the x and y directions 
respectively. We denote the eigenvalues of i7(k) as e„(k) 
and the corresponding eigenvectors as Xn(k, r, s). (Here 
s = 1,2 correspond respectively to fiuctuations of ( and 
If.) Since H is Hermitian and represents a stable ground 
state, the eigenvalues are all positive semidefinite, and 
the eigenvectors are normalized in the usual way, 



X„(k))'x«'(k) = (5„ 



(24) 



We can express a generic fluctuation ip of the phase field 
(j){r) from the ground state at any r inside the unit cell in 
the presence of a unit amplitude of the lowest eigenvector, 
which we call xo(k), as 



ip{r) 



where Nr 



1 



^^e"^"-p„(k)x„(k,r,2), 



(25) 



Lj. Ly 



is the number of unit cells. The reality 
of f forces p* (k) =p„(— k) with the convention that 
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FIG. 9: Spin stiffness as a function of potential strength V, 
normalized to nearest-neiglibor ferromagnetic excfiange J. 



X*(k) = xn(-k). 
easily seen to be 



The energy of such a fluctuation is 



C/ = ^e„(k)K(k)| = 



(26) 



To make the calculation convenient, we choose k^j^ — 
with rain = # being the minimum nonzero 



.min^x with kx 

value allowed. Then, bearing in mind the reality of if, we 
can express this as 



(p(r) 



:Po(kmm)xo(kmm, r, 2) sin I 



(27) 



Comparing the energy of Eq. ((26)) to that of the contin- 
uum version Eq. ()23p , we obtain 



£o(kmiTi) 



(28) 



To be precise, this is the spin stiffness computed at T = 0, 
because thermal fluctuations have not been included. 

The results for the spin stiffness (Fig. |9|) as a func- 
tion of potential strength V show prononounced jumps. 
The discontinuities in JCs coincide with transitions in 
the topological density of the ground-state configuration. 
Notice that the spin stiffness is generally of the order 
of J, the nearest-neigbor ferromagnetic coupling in Eq. 
|3l Depending on the topological density in the ground- 
state configuration, JCs/ J can be either significantly en- 
hanced (for example V/J< 8.15) or reduced (for exam- 
ple, 3.75 < V/J < 5.5). 



IV. FLUCTUATION EFFECTS 

In a purely quadratic theory, fluctuations will have 
little effect on most quantities of interest. However, the 
Hamiltonian of our system contains non-linear couplings 



which can have important qualitative effects. Already 
included in our groundstate analysis is the underlying 
0(3) nature of the spins, which supports the (charged) 
morons that are induced by the periodic potential. Ther- 
mal fluctuations in which meron-antimeron pairs are gen- 
erated above the groundstate can spoil the spin stiffness 
of the system, rendering it dissipative, above an effec- 
tive Kosterlitz-Thouless (KT) transition temperature. In 
principle, interlayer tunneling may have even more pro- 
found effects, for example removing the possibility of 
a true thermodynamic KT transition in the clean limit 
[2^ . In this section we describe how fluctuation effects 
in this model can greatly suppress the impact of inter- 
layer tunneling, and furthermore lead to a lowering of the 
KT transition temperature, particularly near transitions 
between different groundstates. 



A. Suppression of Interlayer Tunneling by 
Collective Modes 

We begin with an evaluation of the effective {i.e., 
rcnormalizcd) interlayer tunneling amplitude due to 
quadratic fluctuations around the ground state. As dis- 
cussed in the Introduction, this can be greatly depressed 
by the presence of the Q-mode, particularly if the gap 
is quite small, as sometimes happens in the vicinity of a 
groundstate transition. 

The fluctuations around the ground state can be ex- 
panded as 

w{r, s,t)^-^Yl ^^A(k)e'(''"— VA(k, r, s), (29) 



where is the number of unit cells in the lattice and 
WA(k) is the amplitude of w in the mode A, k. The reality 
of w means that (V'A(k))* is an eigenvector of iJ(— k) with 
eigenvalue —uJx(h), that is 



V'A(k,r,s 

('^'A(k)) 



■0_A(-k,r,s) 
tZ)_A(-k). 



(30) 
(31) 



Now we are ready to re-express the action in terms of 
these modes. In order to carry out thermal averages, we 
use the imaginary time path integral 

-S^ Jdtfj:^ §[C(r)^(r) - ^(r)C(r)] (32) 



-5 Err',..' ^(r, s)H,s' (r, r')w{r',s')j . (33) 

When we expand in normal modes, w now acquires a 
dependence on iujn as well: 

wir,t)= ^ ^ ziA(k,iw„)e^('^-"-""*)i^A(k,r)(34) 
(wA(k,iw„))* = w_A(-k, -ia;„). (35) 
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The last condition means that we can take wxyoihjiuJn) 
over the entire BZ and for all ia;„ to be independent, 
while the wa<o are dependent. Substituting the expan- 
sion for Sw in terms of the normal modes, and using the 
properties of the wavefunctions Eqs. (|21l22l3ip we obtain 



1 



iw„,k,A>0 



(36) 



As can be seen, the real part of —S is negative definite, 
as it should be for the convergence of the path integral. 
This simple form for the action leads to 



IUJ„ 



l'^A(k) 



(37) 



Now we are ready to find the correlation functions of the 
fluctuations w. The simplest is < (^(r))^ >. Clearly, 
this will be 



1 



{iwiv,s)y) = -— 



iLj„,k,A>0 



|^A(k,r,s)|^ 



-WA(k) 



(38) 



The sum over za;„ can be done by standard methods to 
get n\{]<i) + i, where n\{]^) — (exp/3ajA(k) — 1)"-^ is the 
boson occupation function, so that 



((z«(r,s))2) = ^ |V'A(k,r,s)nn,(k) + i] 



(39) 



k,A>0 



We now consider the effect of these fluctuations on 
inter layer tunneling. The effect as we shall see is most 
pronounced where the Q-mode gap is small. In our pe- 
riodic model this occurs at specific values of V] in a real 
disordered sample, there will always be regions where 
this is true. We add to our original Hamiltonian a small 
intcrlaycr tunnelling term, which may be written in the 
form 



Hilt = -hY,^/^ - '4^? cos^{t). 



We approximate the Q-mode dispersion as 



WQ(k) = E, 



QO 



(40) 



(41) 



We now consider fluctuations of 0(r) which consist of two 
low energy parts, one controlled by the G-mode and the 
other controlled by the Q-mode, under the assumption 
that T is much smaller than the energies of the other 
modes. Wc thus write 



(r) 



0Q(r). 



(42) 



We then integrate out the Q-mode to obtain a "renor- 
malized" interlayer tunneling term. 



HiLT,R = -/^E - ^-^Y cos0G(r)e-<*«W'>/2. 



(43) 



Using Eqs. (|39l) and (|4T|) . we rewrite the last exponential 
as 



exp 



BZ 



^|V'Q(k,r,2)p(nQ(k)-fi). (44) 



For T ^ Eqq one can approximate the occupation num- 
ber of the Q-mode as 



Eqo + ak2 ' 



(45) 



so that as Eqq — > 0, one obtains a logarithmic divergence 
in the k integral, provided ^/'Q(k, r, 2) remains nonzero as 
k — > 0. We have verified that indeed it does, leading to 
a renormalization of h of the form 



ho M he 



(46) 



Thus, in situations where Eqq —s- 0, ft- is strongly sup- 
pressed. 



B. Suppression of Tkt 

Recall that at a generic transition in our periodic po- 
tential model Eqq does not vanish. However, as we will 
now show, the Berezinskii-Kosterlitz-Thouless transition 
temperature Tkt separating the low-temperature, quan- 
tum Hall ferromagnet phase from the high-temperature, 
paramagnetic phase is nevertheless strongly suppressed 
near transitions. Physically, this is because there are two 
local minima with nearly the same energy, but different 
topological density. This makes it very easy for the sys- 
tem to screen vortex/antivortex charge, leading to a very 
small core Ec energy for vortices. 

To see this quantitatively, we introduce and analyze an 
effective, two-dimensional Coulomb gas description of the 
bilayer quantum Hall ferromagnet. At finite temperature 
T, we approximate the partition function of the quantum 
Hall ferromagnet Z « I]{m(R} 6xp(-i?[{m(R}]/fcBT) 
as a sum over the positions of vortices {to(R} with a 
Coulomb gas energy functional 

E[{m{R}]= \ E /C,m(R)m(R')log(^^) 



R^R' 



R 



(47) 



where /C^ is the bare spin-stiffness computed in Section 
IIIIB| m(R) is the vorticity of a plaquette at position R, 
^ is an ultraviolet cutoff scale, and E^ is the core energy 
of each vortex/antivortex, whose value depends on V ^ 
the potential strength, through the potential-dependent 
details of meron core structure. 

To estimate E^iy^, we compare the T = ground- 
state energy (computed in Section [H]) of the two compet- 
ing states whose energies cross at the critical potential 
strength Vc. In the simplest case State 1 is a topologi- 
cally trivial configuration with no vortices, and State 2 
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FIG. 10: Finite temperature phase diagram of the bilayer 
quantum Hall ferromagnet in a periodic potential V (Mea- 
sured in units of J the nearest-neighbor ferromagnetic ex- 
change energy, using a 161b x 16£b unit cell with local 
Coulomb repulsion energy U = 8J). Dramatic decreases in 
the transition temperature Tkt (measured in units of J/ks) 
separating the low-temperature ferromagnetic phase from 
high-temperature paramagnetic phase occur due to changes 
in the topological density of the ground-state at critical po- 
tential strengths. 



has two -|- and two — vortices in a checkerboard pat- 
tern in the unit cell. (See, for example, the configuration 
shown in Fig. [T]) . Generically, they cross with some slope, 
with the energy difference per site being 



AE ^ E2 - El ^ a(y - Vc), 



(48) 



where Vc is the transition point. Wc obtain a « 0.12 
numerically by comparing the energies of the ground and 
metastable states near Vc- Without loss of generality, we 
choose the vortex core size ^ in (|Tf)) such that the vortex 
core energy vanishes as one approaches the transition. 
Together with this choice implies the following form 
for the vortex core energy in the vicinity of the transition 
at Vc-. 



Ec^^^W- 



Vc\ 



(49) 



Note that as one passes through the transition, the high 
and low energy states interchange roles, and (|T7)) still 
governs the fiuctuations of the system for V < Vc, with 
the "absence of a vortex" playing the role of vortices in 
the free-energy functional. Thus the vortex core energy 
Ec{V) is non-negative across the transition at Vc- 

Following the standard analysis @ of the two- 
dimensional Coulomb gas (|47|) . we find an implicit equa- 
tion for Tkt in terms of the vortex core energy Ec and 
the ground-state spin stiffness K-s, 



Ec 



B-l^KT 







In 


\ ' ( 



1 - 



2kBTKT 



(50) 



We plot in Fig. (|10p the result of numerically solv- 
ing (Uni and ([50]) for Tkt as a function of the potential 



strength V. This is one of the central results of this pa- 
per. Note that typically, Tkt is suppressed by an order 
of magnitude compared to the Hartree-Fock predictions 
for a bilayer quantum Hall ferromagnet with no periodic 
potential. 

Near the ground-state transitions, the sharp reduc- 
tion in Ec, the core energy for vortex excitations, de- 
scribed by p9)) . is the principal cause of the dramatic 
suppression of the low-temperature quantum Hall ferro- 
magnet phase. Far from the ground-state transitions, a 
secondary effect of the periodic potential on the phase di- 
agram arises through the ground-state spin stiffness ICg'- 
ksTKT 'r:lCs/2 for \V - Vc\ » inKs/a. As shown in 
Fig. [SI the periodic potential reduces Kg dramatically for 
some values of the potential strength V , leading to a sup- 
pression of Tkt compared to the transition temperature 
in the absence of a periodic potential. 

The Coulomb gas description (|47|) is expected to over- 
estimate Tkt since it contains only a single type of vor- 
tex, and focuses on the spin configurations (States 1 and 
2) with the lowest T = energy. In the critical regime 

ss Vc we observe numerically that the T = energy of 
other states approach that of States 1 and 2. Describing 
the effects of these other states would require generalizing 
(|T7)) to consider multiple types of vortices. On general 
grounds, one expects that including more types of vor- 
tex fluctuations reduces the finite-temperature stability 
of the bilayer quantum Hall ferromagnet phase, so the re- 
sult reported in Fig. [TU]is expected to be an upper bound 
on Tkt- 



V. CONCLUSIONS, CAVEATS, AND OPEN 
QUESTIONS 

It has been clear for some time that there are qual- 
itative discrepancies between theoretical predictions for 
clean systems and the actual phenomenology oi v = \ 
quantum Hall bilayers at small layer separation. There 
is a wide consensus in the community that these discrep- 
ancies are due to quenched disorder. 

Our goal in this paper is to investigate some of the 
nonperturbative effects of quenched disorder in balanced 
V = 1 quantum Hall bilayers. Our principal premises 
are: (i) A strong periodic potential can mimic some of 
the nonperturbative effects of disorder, and (ii) one needs 
to focus only on the (pseudo)spin physics. While these 
premises (especially (i)) can be debated in the context of 
bilayer quantum Hall systems, there is ample historical 
evidence of the fruitfulness of (i) in the Bose-Hubbard 
model and the problem of the quantum Hall plateau 
transition. Other limitations of this work include the 
neglect of long-range Coulomb interactions between in- 
duced topological densities, the neglect of the renormal- 
ization of the classical ground states we found due to 
fluctuations, the neglect of interactions between the low- 
lying modes, especially the G- and Q-modes, and the 
perturbative treatment of the interlayer tunneling. We 
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see this work as a first step in a systematic investigation 
of these effects. 

Starting from the (rather minimal) model Hamilto- 
nian, Eq. we study the ground states and the col- 
lective excitations numerically for interlayer tunneling 
h = 0. As the strength of the periodic potential in- 
creases, we observe generically first-order transitions be- 
tween states with different topological densities. Occa- 
sionally, these transitions are weakly first-order, with a 
new, charge-carrying, quadratically dispersing mode be- 
coming nearly gapless at the transition. Such a mode can 
suppress the interlayer tunneling amplitude strongly at 
nonzero T. Even when the transitions are strongly first- 
order, we show that vortices become very easy to create, 
and drive the Berezinskii-Kosterlitz-Thouless transition 
temperature to zero at the transition. 

ft is important to note that there is a qualitative differ- 
ence between weak and strong h in our approach. While 
our numerics are carried out only for h ~ 0, the inclusion 
of an interlayer tunneling amplitude much smaller than 
any other energy scale will not affect any of the qualita- 
tive physics we uncover. However, if the strength of the 
interlayer tunneling h is increased while other parameters 
such as Ec and the strength of the periodic potential V 
are kept fixed, we expect the system to undergo a set of 
ground state phase transitions leading ultimately to the 
uniform ferromagnetic ground state at very large h. This 
qualitative distinction between weak and strong tunnel- 
ing does not exist in the clean modclQ, but is consistent 
with experiments. 

Based on our results for the periodic potential, we 
can speculate about the effects of true quenched dis- 
order. There are several important effects: Firstly, it 
can be shown via a mapping to the random field Ising 
modcljl^ that disorder converts the sequence of first- 
order transitions we found as V increases into a sequence 
of second-order transitions. While this mapping is rigor- 
ously provable only for classical models, it is believed to 
hold for many quantum phase transitions as well [53, 
(the exceptions seem to be transitions that are com- 
pletely smeared out and destroyed by disorder [s^). This 
sequence of second-order transitions maintains the qual- 
itative distinction between weak and strong tunneling 
mentioned earlier. Furthermore, at the quantum level, 
this implies that the Q-mode generically becomes gapless 
at ground state transitions with true disorder. Secondly, 
even away from a ground-state transition, in a system 
with quenched disorder there will be large regions which 
are close to a transition (the system is in a Griffiths 
phase [iOl). Being in a Griffiths phase also means that 
excitations of arbitrarily low energy are available from- 
large rare regions close to the transition, which leads to 
divergent low- freq uency susceptibilities throughout the 
Griffiths phase [6l| at T = 0. For T ^ 0, since vortices 
are easy to create in such regions, Tkt is expected to be 
strongly suppressed throughout the sample. This is con- 
sistent with the fact that no low-temperature phase with 
the phenomenology of the BKT power-law phase (includ- 



ing a Josephson-like delta-function peak in the interlayer 
conductance) has yet been observed in experiments. This 
suggests two possibilities: (i) The true ground state of the 
system at T = is ferromagnetically ordered, but current 
experiments have probed only the T > Tkt regime, (ii) 
The true ground state is quantum disordered due to some 
combination of quantum fluctuations /quenched disorder. 

Even in the more conventional possibility (i), there are 
several natural sources for the dissipation seen in experi- 
ments. The dissipation could be due to unbound vortices 
in the hydrodynamic transport regime[62j fko <^ k^T, 
and/or due to the low-energy Griffiths Q-modes. 

A second aspect of disorder is that as one crosses do- 
mains of different topological density, the spin wave ve- 
locity varies sharply. This is expected to lead to chaotic 
reflections of spin waves leading to a diffusive-like behav- 
ior at macroscopic length-scales [63j. 

Yet another aspect concerns the critical counterflow 
velocity, which determines the critical counterflow cur- 
rent. In regions where the Q-mode has a very small 
gap, the critical velocit y w ill also be small, and non- 
linear/dissipativc effects [25| will be visible at very tiny 
counterflow current. Similar physics holds for the critical 
interlayer tunneling current density. Some recent work 
takes the point of view that perhaps the experiments are 
in the T < Tkt regime, but the smallness of the critical 
tunneling current density is a primary mechanism driving 
the observed dissipation p^. 

Let us now turn to issues of zero-temperature physics. 
It is possible fluctuations quantum disorder the ground 
state (even in the periodic-potential model) near the 
mean-fleld ground state transitions that we found. 
Clearly, one needs an effective theory of the low-lying 
modes and their interactions to address these issues. 
An important source of quantum disordering is the tun- 
neling of (multiple) vortices. Due to the spin-charge 
relation, the quantum tunneling events (called hedge- 
hogs) of a single vortex will violate charge conserva- 
tion. However, especially near transitions, one can imag- 
ine multiple vortices in a unit cell with total vortic- 
ity zero tunneling together. This has a connection to 
ideas of deconfincd criticalitv [63 - l66| . In this type of sce- 
nario, first described for quantum dimer[63| models and 
quantum antiferromagnets[65| . dimers/spin excitations, 
are described as composites of monomers/spinons. Sin- 
gle hedgehog events are forbidden due to lattice sym- 
metries, but multiple hedgehogs are allowed. In cer- 
tain well-defined extentions of quantum antiferromag- 
nets, the multiple- hedgehog events can be shown to be 
irrelevant [671. leading to a critical point with deconfined 
spinous [66[. 

The important difference is that in the picture of the 
quantum Hall bilayer presented here, it is not lattice sym- 
metries but the spin-charge relationjs^, which en- 
forces the absence of single hedgehog events. Thus, if 
deconfined criticality were to occur for the periodic po- 
tential model, it would also likely occur in the model with 
true disorder. 
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The quantum Hall bilayer remains a rich system which 
potentially supports a host of physical behaviors yet to 
be explored, particularly in the presence of strong dis- 
order potentials that are almost certainly a feature of 
their realization in semiconductor systems. We believe 
immersing the system in a periodic potential offers a win- 
dow through which one may begin an exploration of this 
physics. 
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